#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Thu Jun 20 18:53:27 2024

@author: ndrach
"""

import numpy as np
import matplotlib.pyplot as plt

V1, V2, V3 = 150, 400, 1000
rc = np.logspace(-9,-6,100)
d = 5e-3


gamma = 0.073
e0 = 8.85e-12
ER = 2*np.sqrt(gamma/(e0*rc))
ER2 = np.sqrt(np.pi*gamma/(e0*rc))
ER3 = 1/np.pi * np.sqrt(3*gamma/(e0*rc))



fig, ax = plt.subplots(figsize = [3.5,2.75])

ax.semilogx(rc,ER*1e-9, lw = 2, color = 'k', label = 'Rayleigh limit')

ax.plot(rc,np.ones(len(rc))*0.75, ':', color = 'k')#, label = 'ion evap. threshold')

ax.set_ylim(0,2)
ax.set_xlim(3e-9,1e-6)
ax.legend(fontsize = 6)
ax.set_ylabel('Electric field strength (V/nm)', fontsize = 8)
ax.set_xlabel('Tip radius (m)', fontsize = 8)
ax.tick_params(labelsize = 8, which = 'both', direction = 'in')

